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ABSTRACT 



Colliding winds of massive star binary systems are considered as potential 
sites of non-thermal high-energy photon production. This is motivated merely 
by the detection of synchrotron radio emission from the expected colliding wind 
location. Here we investigate the properties of high-energy photon production 
| in colliding winds of long-period WR+OB-systems. We found that in the domi- 

nating leptonic radiation process anisotropy and Klein-Nishina effects will likely 
yield spectral and variability signatures in the 7-ray domain at or above the 
sensitivity of current or upcoming 7-ray telescopes. Individual colliding wind 
binary (CWB) systems are therefore regarded as candidate 7-ray sources for the 
satellite-based 7-ray missions AGILE and GLAST, and ground-based imaging 
atmospheric Cherenkov telescopes (MAGIC) or telescope arrays (CANGAROO, 
H.E.S.S., and VERITAS). 

Analytical formulae for the steady-state proton- and electron particle spectra 
are derived assuming diffusive particle acceleration out of a pool of thermal wind 
particles, and taking into account adiabatic and all relevant radiative losses (In- 
verse Compton scattering, synchrotron radiation, non-thermal bremsstrahlung, 
Coulomb losses and inelastic proton-proton collisions). For the first time we 
include their advection/convection in the wind collision zone, and distinguish 
two regions within this extended region: the acceleration region where spatial 
diffusion is superior to convective/advective motion, and the convection region 
defined by the convection time shorter than the diffusion time scale. We show 
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that electron losses my well extend into the Klein-Nishina transition regime, and 
find analytical approximations for the loss rate in this regime. 

The calculation of the Inverse Compton radiation uses the full Klein-Nishina 
cross section, and takes into account the anisotropic nature of the scattering 
process. This leads to orbital flux variations by up to several orders of magnitude 
which may, however, be blurred by the geometry of the system (eccentricity, 
inclination). The maximum inverse Compton radiation occurs when the WR- 
star is located behind the OB-star. Both, anisotropy and Klein-Nishina effects 
will likely yield characteristic spectral and variability signatures in the 7-ray 
domain. Non-thermal bremsstrahlung emission is found to be mostly of minor 
relevance. Propagation effects lead to a deficit of low-energy particles in the 
convection zone, which - if its size is sufficiently large - may leave visible imprints 
in the radiation spectra. If protons are accelerated to at least several GeV, 7r°- 
decay 7-rays contribute to the high energy SED, and, depending on the injected 
electron-to-proton ratio, might be visible with upcoming 7-ray instruments. We 
show that pair production from photon-photon collisions can not be neglected in 
these systems in general, and might affect the emitted spectrum above ~ 50GeV 
depending on orbital phase and system inclination. 

The calculations are applied to the typical WR+OB-systems WR 140 and WR 
147 to yield predictions of their expected spectral and temporal characteristica 
and to evaluate chances to detect high-energy emission with the current and 
upcoming 7-ray experiments. 

Subject headings: Stars: early-type - Stars: binaries - Stars: winds, outflows - 
Gamma rays: theory - Radiation mechanisms: non-thermal 

1. Introduction 

Early type stars (0-, early B-, Wolf-Rayet (WR) stars) are hot stars (T cfr > 10000 K) 
with masses > 20 M Q . They are known to possess some of the strongest sustained winds 
among galactic objects. As a class Wolf-Rayet stars have the highest known mass loss 
M ~ 10~ 4 -~ 5 M /yr of any stellar type. The supersonic winds of massive stars may 
reach terminal velocities of > 1000 — 5000 km/s (Cassinelli 1979), their kinetic energy 
L w = however, rarely exceed 1% of the bolometric radiative energy output of 

typically ~ 10 38 erg/s. 

In recent years massive stars have been connected to several high energy phenom- 
ena: they are suspected to be the progenitor of some type of 7-ray bursts (Woosley 1993; 
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Paczynski 1998), and are known to be an interacting medium for the blast wave expelled 
from supernova explosions. Triggered by the detections of several 7-ray sources by EGRET 
that are not unambiguously identified yet but found in positional coincident with massive 
binary systems (Kaul & Mitra 1997; Romero et al. 1999), and motivated by the detection 
of synchrotron radiation from the collision region in some massive binaries, these systems 
have also been proposed as potential sites of non-thermal high-energy photon production 
(Montmerle 1979; Eichler & Usov 1993; Romero et al. 1999). 

Unlike in a wind of a single massive star where particles has been proposed to be ac- 
celerated by either multiple weak shocks from line-driven instabilities (Lucy & White 1980; 
White 1985) to compensate for the expansion and radiative losses even close to the stellar 
photosphere, or at the terminal shocks created by the interaction with the swept up inter- 
stellar medium (Casse & Paul 1980; Volk & Forman 1982), the collision of supersonic winds 
produces strong shocks where both electrons and protons can be efficiently accelerated to 
high energies through first order Fermi acceleration (Eichler & Usov 1993). 

In this latter scenario the shock region is exposed to both, a strong radiation field 
in the UV range from the participating hot stars, and their magnetic fields. The detec- 
tion of synchrotron radiation from such collision regions, implied by the flat or negative 
spectral indices and brightness temperatures of 10 6 " 7 K, proofs the existence of magnetic 
fields as well as relativistic electrons in the collision region of some massive binary systems 
(Abbott et al. 1986). In addition, observations of radio emission reveal an extended region 
of the synchrotron radiation (Dougherty et al. 2000). In fact, in some binaries (e.g. WR 147, 
WR 146, OB2 No. 5) the collision region between the main sequence stars has been resolved 
in the radio band showing an extended, slightly elongated non-thermal feature on VLA and 
MERLIN images (Dougherty et al. 1996; Dougherty et al. 2000; Contreras et al. 1997) in 
addition to the free-free emission from the spherical wind of the stars. Recently, the ex- 
tended wind-wind collision region from WR 147 has even been detected in the X-ray band 
by Chandra (Pittard et al. 2002). 

Those electrons with energy r y e m e c 2 will inevitably also be responsible for a non-thermal 
high energy component at ~ l1 er v produced through inverse Compton scattering of the 
dense stellar radiation field with characteristic energy ex- Interestingly, recent XMM and 
simultanous VLA observations of the WR-star 9 Sgr might already lurk the hint of a non- 
thermal X-ray component. While the hard X-ray component of 9 Sgr could be equally well 
fitted by either a hot multi-temperature thermal model with kT > 1.5 keV, or a steep power- 
law with power-law index > 2.9, suggesting a compression ratio < 1.8, the VLA spectrum 
at the time of the XMM-observations was clearly non-thermal, thus indicating a similar 
compression ratio ~ 1.7 (Rauw et al. 2002). 
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The relativistic electrons will also loose some fraction of their energy due to non-thermal 
bremsstrahlung in the field of ions that are embedded in the stellar winds. It has been 
estimated, however, that the resulting non-thermal bremsstrahlung component in the 7-ray 
domain at ~ ( , y e /2)m e c 2 from electrons of energy r y e m e c 2 will be a minor contribution to the 
overall emission because of the ambient gas densities involved (Benaglia & Romero 2003). 

If electron acceleration is also accompagnied with acceleration of protons out of the 
thermal pool of wind material, then proton interactions with the ambient ions will produce 
7-rays through the hadronic neutral pion decay channel p+p — > tt°+X, 7r° — > 7+7- Addi- 
tionally, radiation from the secondary pairs, generated through the decay of charged mesons 
that are produced by the hadronic pp-collisions, is expected to contribute also to the overall 
broad band spectrum. The 7r°-decay radiative channel has already been considered in the 
past in the context of winds from massive stars ((Chen & White 1991; White & Chen 1992); 
White 1985, for the case of single O-stars), however, no propagation effects inside the ex- 
tended collision region nor competing loss mechanisms (such as expansion losses in the wind) 
has been taken into account here. 

The goal of this work is to extend the model of non-thermal emission in the high-energy 
domain from the colliding wind region of binary systems of massive stars to include all 
relevant energy losses and simultanously also consider the propagation effects that affect 
a relativistic particle distibution in such an environment. After a brief discription of the 
geometrical model considered in this work (Sect. 2), we evaluate the evolution of both, 
proton and primary electron spectra, on the basis of a simplified diffusion-loss equation. In 
Sect. 3 we calculate the expected photon emission due to the inverse Compton process (for 
the first time including Klein- Nishina and anisotropy effects) , relativistic bremsstrahlung and 
the 7-rays from the decay of n° that is produced in hadronic proton-proton collisions. The 
total expected gamma-ray spectrum is corrected for photon absorption in the UV radiation 
field of the massive stars. We apply our model to the famous long-period binary system 
WR 140 and the above mentioned WR 147 in Sect. 4. Our conclusions summarize our results 
and provide an outlook on the detectability of colliding wind binary systems with upcoming 
instruments like GLAST, and contemporary Imaging Atmospheric Cherenkov Telescopes 
(IACTs). 



2. The geometric model of a colliding wind region 

The typical radial velocity profile V(x) of winds from hot massive stars obeys the relation 

V(x) = v 00 {l-r a /xf 
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(Lamers & Casinelli 1999) with (3 ~ 0.8, r s the radius of the star and Voo the termi- 
nal velocity. Recent observations indicate the existence of clumping in the wind (e.g. 
Moffat 1996; Schild et al. 2004). For our schematic picture here we shall postpone the ef- 
fects of clumping to later work and will consider a homogeneous uniform wind. We also 
neglect any small-scale shocks within the wind. The winds from binary systems flow nearly 
radially until they collide to form a discontinuity at the location of ram pressure balance, 
and forward and reverse shocks. In the shock region the gas is heated to temperatures of 
10 6 - 8 K (Luo et al. 1990; Stevens et al. 1992; Usov 1992) which causes strong thermal X-ray 
emission. Behind the shock the gas expands sideways from the wind collision region along 
the contact surface out to larger r. It is therefore suggestive to distinguish two regions of 
the extended emission site (see Fig. 1): In the acceleration zone r < ro (first-order) diffusive 
acceleration provides high-energy particles out of a thermal pool. While the stellar winds pro- 
hibit the escape of energetic particles on the upstream side of the shocks, we anticipate that 
in the downstream region spatial diffusion is more efficient than convective/advective motion 
(which we call "convection" in the following) in transporting particles to the boundary of 
the acceleration zone at r = r . Subsequent to their diffusive escape from the acceleration 
zone, the energetic particles enter the convection dominated zone r > r . The characteristic 
radius r = r is defined by equality of the convection time scale t con v and the diffusion time 
scale tdifr- 

Furthermore, we demand that the distance of the emission region from the low- 
momentum-wind source is large compared to the longitudinal extension of the emission 
region. While in reality the form of the contact surface is bent towards the star that 
shows the lower wind momentum, observed as arcs of emission on radio images (e.g. 
Dougherty et al. 2005), we consider here the simplified geometry of a cylinder disk (see 
Fig. 1). This may be justified by the rapid particle energy loss rates that do not allow 
the transport of energetic particles to large distances from their acceleration site. In other 
words, most particle energy losses are expected to appear close to the acceleration zone 
where a cylinder geometry appears a reasonable approximation. The thickness d of the 
cylinder-like emission region is governed by the velocity of the hitting winds. We shall 
further neglect here the interaction of the stellar radiation fields on the wind structure 
(Gayley et al. 1997; Stevens et al. 1994) which is justified for long-period binaries. 

In the case of a collision between the spherical wind of a primary (e.g. a WR-star) and 
a secondary companion (e.g. a OB-star) which has reached their terminal velocities (vob, 
the location of the shock is determined by the balance of the ram pressure of both 

winds 

Pob^ob = Pwr^wr (!) 
where pqb, Pwr are the densities of the gas ahead of and near the shock of the stellar winds 
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of the OB- and WR-stars. In this case the distance of the shock front from the WR-star, 
Xwr, and from the OB-star, xqb, is given by 

1 /rj 
xwr = — — -=D, x B = — — -=D (2) 

i + Vv i + Vv 

where 

MOB^OB /q\ 

rj = -. (3) 

M WR t> WR 

and D = iwr + ^ob is the separation of the binaries and Mob, Mwr are the mass loss 
rates of the WR- and OB-star, respectively. Since M b < M WR and vob ~ %r, the shock 
location is rather close to the OB-star, i.e. D 3> x = xob- This appears to be in excellent 
agreement with the observed locations of the collision region in the radio domain of e.g. WR 
147 (Dougherty et al. 2000), and this also implies that the shocked gas of the winds of both 
the OB- and the WR-star have a comparable mass density. 

If the wind collision occurs at a distance smaller than the Alfven radius 

rA ~ r s x (1 + for £ < 1, 
r A » r s x for ^ >^ j 

(where £ = B 2 s r 2 s /{2Mv O0 ), B s is the star's surface magnetic field) from the OB- 
star, significant deceleration of the WR-wind in the OB-radiation field is expected 
(Eichler & Usov 1993). In this case, the colliding winds do not reach their terminal ve- 
locities. Typically, £ ~ 0.01... 0.1 leading to ~ 3 . . . 5r s for OB-stars (Barlow 1982). 
For the present work we limit ourselves to situations, where the binary winds reach their 
terminal velocities at the shock location. 

The external magnetic field B of a star changes in an outflow from the classical dipol 

field 

B ss B s x {r s /x) 3 
to a radially dominated one at the Alfven radius 



B S3 B s x r 6 Jr A x 

and finally to a toroidally dominated one for x > r s Voo/v ro t 

B S3 B s x v TOt r 2 J (vooTax) 

where v rot is the surface rotation velocity with typical values of order ~ 0.1i>oo for early 
type stars (Conti & Ebbets 1977; Penny 1996). We use these equations to determine the 
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magnetic field at the location of the collision front. This value is held constant throughout 
the emission zone. 

The surface magnetic field of massive stars are not well known. (Donati et al. 2002) 
report the detection of a 1.1 kG dipolar magnetic field of presumably fossil origin at the 
surface of the young O star 6>i Ori C. (Ignace et al. 1998) argues for surface magnetic 
field strengths of order 10 4 G in WR-stars. On the other hand from the non-detection 
of the Zeeman effect for many O- and B-stars only upper limits of order a few 100 G ex- 
ist (Barker 1986; Mathys 1999). For this work we shall fix the surface magnetic field at a 
reasonable value of B s = 100 G unless stated otherwise. 

The plane of the binary system is inclined by an angle % with respect to the observer 
(i = 7i / 2 corresponds to an observer in the plane of the stars), and <pB, the angle between 
the projected sight line and the line connecting the stars is a measure of the orbital phase of 
the system. In the following, periastron passage is defined by the orbital phase $ = 0, and 
4>b = for the WR-star being in front of the OB-star along the sight line. 

In the co-rotating system centered on the OB-star the location of the emission site is 
defined in polar coordinates by the azimuthal angle and the polar angle 9 = arcsin(r/x) pa 
r/x for r <i (see Fig. 1). In the same star-centered frame the line-of-sight to the observer 
is described by the angles 0l with tan0^ = sin0s coti, and 9l with cos#l = cos 0b sini 

The radiation yield of inverse Compton scattering depends on the angle # p h = 9 SC be- 
tween the directions of the incoming (from the OB-star) and outgoing photons, which ob- 
viously depends on the location of the scattering electron as well as the orbital phase. We 
find 

/U P h = cos 9 ph = cos 9 L cos 9 + sin 9 L sin 9 cos(0 — <p L ) 
The corresponding azimuthal angle p h is irrelevant for the scattering process. 

3. Particle spectra 

We expect two standing shocks and a discontinuity between them. For typical massive 
stars the wind velocities are of comparable value and in addition the ram-pressure balance 
(Eq. 1) ensures a similar upstream gas density for both the OB- and the WR-wind shocks. 
It thus appears justified as assume the two shocks as well as the corresponding acceleration 
and escape rates to be identical. 

As motivated above we distinguish two zones of the extended emission region. In 
the acceleration zone suprathermal particles of energy E from the stellar winds are en- 
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ergized through diffusive shock acceleration at a rate E = aE where a = V, 2 (c r — l)/(3c r K a ) 
(Schlickeiser 2002) with V s = t>oB = %R is the upstream velocity of the standing shock, c r 
the compression ratio and n a the (assumed) energy-independent diffusion coefficient perpen- 
dicular to the wind contact surface. In this region diffusion dominates over the convective 
flow along the wind contact surface. The acceleration region can in good approximation 
be treated as a leaky box with a free escape boundary at r (see Fig. 1). Introducing an 
energy-independent escape time T = rl/(AK d ), where K d is the diffusion coefficient along the 
wind contact surface, the continuity equation for this zone then reads: 

±(EN(E)) + ^ = Q 6(E-Eo) (4) 

where E = aE—Ei oss includes energy gain through diffusive acceleration as well as continuous 
energy losses. The solution N(E) consists of a power law that is modified at the high-energy 
end of the spectrum: N(E) oc E~ s f c (E) where f c {E) depends on the radiative energy losses 
employed, and s = (aT ) _1 + 1 where energy losses are negligible. 

Adjacent to the acceleration zone is the so-called convective zone, where convection along 
the wind contact surface is a faster transport process than is diffusion. For the convective 
flow we assume for simplicity a constant velocity, V = pvoB with p = const < I. The 
continuity equation for this zone also includes adiabatic losses, oc EV V, and is given by 

V(V N(E,r)) + -^[(E- f V V) N(E,r)\ = (5) 

where E = —E\ oss represents the continuous energy losses in this region. The boundary 
conditions N(r — > oo) = and N(r = r , E) = N acc (E) apply, where iV acc (i?) is the 
homogeneous particle density in the acceleration region. 

At the location r = r the diffusive escape time scale T = rl/AK, d equals the convection 
time scale £ conv = r /V, and this allows to determine r = For the power-law index 

s one then finds 

s = 1 + p 2 Sc r ^ (6) 

V 4(c r -I) n d 1 ; 

It is remarkable that for isotropic diffusion, K a = K d , s is fully determined by the compression 
ratio c r of the shock and the ratio of the convection velocity to the shock velocity. For strong 
shocks, c r = 4, hard particle spectra with s = 1+p 2 are then expected in the acceleration zone 
while the generic s = 2-spectrum requires non- isotropic and/or energy-dependent diffusion, 
or an extreme value for the convection velocity V = v Q b- The smallest possible size r of the 
acceleration region corresponds to the diffusion coefficient K d approaching the Bohm-limit. 
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3.1. Electron spectra 

The general solution of Eq. 4 is given by 

dE" 



N(E) = l [ dE'Q(E')exp -1 f 
E J Jo Je> 



E(E") 



(7) 



It has been shown that inverse Compton scattering, if treated in the Thomson regime, 
in most cases determines the maximum electron energy and is the most important ra- 
diative loss channel for ultrarelativistic electrons in colliding winds of massive stars 
(Eichler & Usov 1993; Miicke & Pohl 2002), since typically Wph/r/uB > 1 (w P h,T is the energy 
density of the stellar radiation field, ub the magnetic field energy density in the emission 
region). For typical system parameters one finds w p h,T/ M B ~ 67L bo i 38 / BqXq B13 where Lboi,38 
is the bolometric luminosity of the OB-star in 10 38 erg/s, Bq is the magnetic field in the 
collision region in Gauss and xob,i3 = ^ob/(10 13 cm). At lower energies bremsstrahlung 
and Coulomb scattering determine the shape of the electron spectrum. Radiative losses 
included in our calculations for the electron spectra are synchrotron losses, inverse Comp- 
ton losses on the stellar radiation field of differential photon density n p h,T(e) = no5(e — ex) 
(mono-chromatic approximation), electron-ion bremsstrahlung and Coulomb losses. The 
total radiative energy loss rate is then given by 



6l C,TL = ,-^3 OT U ph 



E = -(b syn + &i C ,tl)£ 2 + (a - b hr )E - b conl (8) 

with 

for the Thomson regime with c the velocity of light and Ot = 6.65 • 10~ 25 cm 2 the Thomson 
cross section, 

h 4 
2 

^br = — Ot Ot C N h 
IX 

in the weak-shielding limit where we have neglected the logarithmic term, and where a is 
the fine structure constant, m e the electron mass and Nh the thermal ion density (in cm -3 ), 
and (Schlickeiser 2002) 

&coui = 55.725 c ot N h m e c 2 . 



Fig. 2 shows the energy loss rates in comparison to the acceleration rate for a set of 
parameters typically found for colliding wind binaries. Obviously Coulomb losses dominate 
the low energy end of the electron spectrum, and simultanously provide a lower limit to 
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the acceleration rate a. At energies E > mlc^/er Klein-Nishina effects significantly modify 
the energy loss rate due to inverse Compton scattering. Fig. 2 demonstrates that firstly, 
Klein-Nishina effects can in general not be neglected in the environment of colliding massive 
star winds. Secondly the curvature of the loss rate due to the Klein-Nishina decline of the 
inverse Compton cross section actually starts already at much lower electron energies. Con- 
sequently, in some cases the electron spectrum might be rather limited by synchrotron losses 
if the wind parameters are favourable. In the following, we take into account Klein-Nishina 
effects already above E > E TL where E TL is considered the energy below which the Klein- 
Nishina energy loss deviates not more than 20-25% from the Thomson limit approximation. 
Typically, for the systems considered in this work _E TL 10~ 2 MeV/e T ,MeV with e T ,MeV the 
target photon energy in MeV. 

In the extreme Klein-Nishina regime the total electron energy losses in the acceleration 
zone at high energies are effectively only provided by synchrotron radiation. In the transition 
range between the Thomson and Klein-Nishina limit, however, inverse Compton scattering 
may still dominate over synchrotron emission, thus mandating a proper treatment of the 
former. To allow for analytical solutions of Eq. 4 we approximate the Klein-Nishina decline 
of the loss rate in the transition range by 

£ic,kn = -6ic,TL^ 2 (l-^e T ,Mev/5-10- 2 MeV) (9) 

for all energies E < E s = 10 _1 - 7 - _1 - 6 MeV/e T ,MeV, and 

E ICiK N ^ -(0.27<7 T cn ph!l 4^ eV £ + 7e^ 1 c y 46 * + 2.1 • 10" 4 ) MeV/s (10) 

= -(q a E + q b ) (11) 

for E > E s and with n p h,T the integrated target photon density and D14 the binary separa- 
tion in 10 14 cm. These approximations are tested for ex = 1 • • • lOOeV, and are suitable for 
0.01-u p h,T < «b < M P h,T (while they deviate by more than an order of magnitude at energy 
q a /(2b syn ) [(1 + a/(1 - 46 syn g 6 /g2)] if u B decreases to < 0.003w phjT ). 

This approximation is indicated in Fig. 2 by the long-dashed line and is suitable for 
cases where relativistic bremsstrahlung losses are much smaller than synchrotron or inverse 
Compton losses. While this approach takes reasonably into account the early deviation of 
the Klein-Nishina cross section from the Thomson cross section, it overestimates the energy 
loss rate in the extreme Klein-Nishina regime, which is mainly responsible for the steepness 
of the declining electron spectrum. Consequently, our derived electron and photon spectra 
shall be regarded with caution in their declining part. 

The analytical solutions of Eq. 4 for electrons in the acceleration region are detailed in 
App. A. Fig. 3 shows examples of the resulting electron spectra for various values of the 
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binary separation, while all other parameters are the same as used in Fig. 2. The general 
shape is determined by the interplay between acceleration gains and losses. At low energies 
Coulomb losses dominate and lead to the well-known upturn towards low energies due to 
the stronger increase of the acceleration rate with energy E acc oc E when comparing to 
the competing radiative losses E\ oss oc const. For a reasonable convection velocity V = 
l/2t>oB ( e -g- Luo et al. 1990) and assuming strong shocks c r = 4, a oc E~ 2 particle spectrum 
develops in the acceleration region if n a / K d = 4, which we shall use in the following for 
demonstrative purposes if not noted otherwise. Thus in Fig. 3 a _E~ 2 -spectrum develops 
until inverse Compton losses cause a steepening. The shape of this decline reflects the 
approximation employed to simulate the losses in the Klein-Nishina regime. The kink at 
~ 10 3 5 MeV corresponds to the transition from the Thomson regime to the Klein-Nishina 
loss rate approximation. Finally, the cutoff at energy E c can be either due to radiative losses 
dominating over the acceleration rate, or the diffusion coefficient K d approaching the Bohm 
limit /tBohm ~ l/3r^c where tl is the Larmor radius. In the latter case the escape time scale 
T decreases with energy. For simplicity (and in order to keep the solutions analytical) we set 
T = at K d < /tBohm which then causes a sharp cutoff. More sophisticated calculations may 
smooth this decline towards an approximately exponential shape (Protheroe & Stanev 2000). 
In the former case the cutoff is determined by the synchrotron loss rate that dominates 
regarding the flattening loss curve at high energies. In Fig. 3 the particle spectra cut off due to 
radiative losses for D < 10 14 cm while at larger binary separations the Bohm limit causes the 
cutoff. Here both, synchrotron and inverse Compton losses, cease to be able to significantly 
affect the acceleration spectrum. For comparison we have also calculated electron spectra 
assuming isotropic diffusion. Fig. 4 shows the resulting spectra where all parameters are 
unchanged with respect to Fig. 3 except for the ratio of the diffusion coefficients K a / K d . As 
expected hard particle spectra with spectral index s = 1 + p 2 = 1.25 develop, and radiative 
losses alter the spectral shape depending on the binary separation. For E — > E c the spectra 
diverge for ((a - b hl ) * T ) _1 < 1 (see Eq. A2, A4), which is fullfilled here. 

Fig. 5 shows the energy loss rates in the convection zone for the same parameters as 
used for Fig. 2. On account of expansion in the cylindrical convection flow the magnetic 
field strength and the gas density fall off quickly with distance from the acceleration region. 
Coulomb interactions, bremsstrahlung, and synchrotron radiation will therefore loose impor- 
tance as energy loss mechanisms and the inverse Compton scattering effectively provides the 
bulk of the electron energy losses. 

Using a constant convection velocity the steady-state particle spectrum in the convection 
region can be found by solving Eq. 5. With the boundary conditions N(E, r) — > for r — > oo 
and N(E, r) — > N acc (E) for r — > r where N acc (E) is the steady-state particle spectrum in the 
acceleration region we derive the following analytical solutions using the standard method 
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of characteristics: For E < _E T l the radiative losses are E = — 6 syn &ic E 2 and 



N(E) 



roy/3 



1 | 3& synfcIC £r 



2V 



M 2/3 _i 

r / 



n -2 



iV acc (^) 



with 



E = E 



1/3 



1 + ^^^ 



r ) 



2/3 



- 1 



(12) 



(13) 



For E > E^l an analytical solution can not be found using Eq. 9. In the convection region the 
approximation for the Inverse Compton energy losses in the Klein-Nishina regime, E > Etl, 



Eic,kn = —q c E 
for E TL < E < E cxtr with q c m bi C E TL , and 

^IC.KN = —Qd 



(14) 



(15) 



for E > -Ecxtr = Qd/Qc with qd 9/16caTmlc 4: n p h,Te T 1 , turns out reasonable if -Etl is 
increased to Etl ~ 10 _1 ' 7 MeV/eT,MeV (see Fig. 5). In this case the solution of Eq. 5 is 



iV(E)=exp[|(r-ro)] {^) 2 " iV acc (E) 



with 



£ = £exp[^(r-r 



)] (- 



1/3 



for E TL < E < E cxt , and 



with 



for E > E cxt . 



E = E 



N(E) = (^) N acc (E) 
3q, 



1/3 



+ 



41/ 



1/3 



r - r 



(16) 

(17) 
(18) 
(19) 



The total solution is then combined such to assure continuity for all functions E(r), 
N{E, r) r Q) N &cc {E)). 

Fig. 6 shows the resulting electron spectrum in the convection region at r = r . . . 10 14 cm 
with a binary separation of D — 10 14 cm and using the parameters as described in Fig. 2 
for the acceleration region. As the particles convect along r radiative losses alter the high 
energy end of the particle spectrum leading to a decreasing cutoff energy with increasing r 
while adiabatic losses lower the overall particle density. In Fig. 7 the binary separation has 
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been enlarged to D — 10 15 cm. As a result radiative losses show a significant impact on the 
spectral shape only at large r. 

In summary, taking into account convection in the extended emission region alters the 
energy cutoff in the integrated particle spectrum if radiative losses prove to be important, 
and simultanously lowers the total non-thermal particle density. This behaviour is reflected 
in the corresponding photon spectra (see Sect. 4). 



3.2. Proton spectra 

If protons are accelerated together with the electrons, hadronic nucleon-nucleon inter- 
action take place with an approximated energy loss rate of 

E = -b pp E (20) 

above the kinematic threshold for pion production E > E thr ~ 0.28 GeV 
(Mannheim & Schlickeiser 1994), 

b pp = 1.3 x 3cN H a pp m n /m p (21) 

where a pp = 3 • 10~ 26 cm 2 is the corresponding hadronic cross section, and the factor 1.3 takes 
into account the here assumed metallicity (90% H, 10% He). Note that this linear relation 
for the energy loss rate is exact only above a nucleon kinetic energy of a few GeV, whereas 
at lower energies it shall be considered as an approximation. The corresponding error drops 
below ~ 30% above ~ 10 GeV. In addition, Coulomb-losses in the dense wind of the massive 
companion star of the WR-star will alter the injection spectrum. The Coulomb-losses are 
(Mannheim & Schlickeiser 1994) 

■ 3 ca T m e c 2 Z 2 In A Ar f3 2 

E = 2 < 22 > 



with 



and 



,JE (E + 2m p c 2 ) 
E + rripC 2 



x m = 0.2 



T, 



e 



10 s K ' 

The Coulomb-barrier occurs at E m ~ 20MeV(T e /10 8 K) where T e « 10 8 K is the electron 
temperature. We approximate the Coulomb-losses by 



E — — 6bcl-EMeV 



(23) 
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below the Coulomb-barrier with 6 be i = 6.6 • 10 16 Z 2 N H /xf n MeV/s, E Mc v is the particle 
kinetic energy in MeV, and 

E = -=^= (24) 

V -&MeV 

above the Coulomb-barrier with 6 a b = 6.7- 10~ 12 Z 2 Nh MeV/s. This approximation deviates 
at most a factor ~ 2 (by approaching the Coulomb-barrier) from the exact loss formula. In 
the relativistic regime, j3 > 0.5, we use 

E = -6 rel (25) 

with b rcl = 3.1 • 10~ 13 Z 2 N H MeV/s. 

Eq. 4 describes the behaviour of the steady-state spectrum in the acceleration region. 
The solution, Eq. 7, can again be derived analytically (see App. B). 

Fig. 8 shows examples of steady-state nucleon spectra with varying distance of the 
binary stars to each other. Close binaries show an upturn at low energies due to a high 
rate of Coulomb losses, and a spectral shape at larger energies that repeats the acceleration 
spectrum due to the same energy dependence of losses from hadronic pp-interactions and 
energy gain. For large binary separations Coulomb losses are unimportant, and the loss 
corrected particle spectrum has the same spectral shape as the steady-state acceleration 
spectrum. Due to the low hadronic cross section radiative losses hardly cause any cutoff in 
the particle spectrum. Instead, faster particle escape when approaching the Bohm diffusion 
limit leads to a steepening of the emitting particle distribution. For simplicity, we chose 
to treat this effect analog to the electron acceleration in Sect. 3.1, i.e. we set Tq = for 

K d ^ K Bohm- 

Applying mass conservation the continuity equation implies Nh oc r~ l for the target 
ion density Nh- With r > tq, where ro indicates the transition from the acceleration to 
convection region with typically r ~ 10 12 cm in the star systems considered here, hadronic 
pp-interactions and Coulomb-losses can readily be neglected. Eq. 5 can be solved to give the 
analytical solution for a constant convection velocity V: 

N(E) = (^) 2/3 N acc (E) (26) 



with 



E E ( - j (27) 



Fig. 9 shows the resulting nucleon spectrum for the same parameter set as used in 
Fig. 6. Obviously the losses due to pp-interactions are of minor importance as compared to 
the adiabatic losses the nucleons suffer in the convection region. 
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3.3. Particle spectra normalization 

For applications to massive binary systems the normalization of the relativistic particle 
component is limited by several constraints: Firstly, the inverse Compton, bremsstrahlung 
and 7r°-decay 7-rays must not overproduce any observational limits imposed by 7-ray obser- 
vations (e.g. EGRET, ...). Secondly, since the particles are accelerated out of the pool of 
thermal particles, particle number conservation dictates that the relativistic particle flux in- 
jected into the system must be smaller than the wind particle flux entering the acceleration 
zone, i.e. Qo < M /(m p x 2 4:n d) with d the thickness of the acceleration site. Further- 
more, due to energy conservation the total injected particle energy can not be larger than 
the total kinetic wind energy of the binary system which rarely exceeds 1% of the total 
radiative energy output of the stars, (typically L w < 10 erg/s). The energy density of 
accelerated particles is given by U acc ~ QoE Q /a assuming an _E~ 2 -spectrum (aT = 1). 
In equilibirium the total acceleration power equals the loss power due to escape, lead- 
ing to L acc = UsjccVsjcc/Tq where V acc = r^nd is the acceleration volume. By noting that 
(ro/4x 2 )L w is the kinetic power available to the acceleration region, Q can not exceed 
Qo ~ (^o/4^ 2 )-^w/(Kcc-E'o) = L w /(4x 2 irdE ). These latter two arguments pose to date a 
stronger constraint on the normalization than the 7-ray limits from EGRET. In the follow- 
ing we use this maximum possible injection power unless stated otherwise. 

4. Photon spectra 

This work is devoted to photon emission at high energies with emphasis on the > 1 MeV 
regime. Important non-thermal continuum emission processes here are inverse Compton 
scattering of the dense stellar radiation field to high energies, relativistic bremsstrahlung of 
the electrons in the field of the ions in the wind and the decay of 7r° — > 7 + 7 that are 
produced in hadronic nucleon-nucleon collisions. For the calculations of the photon spectra 
we assume the particle distributions to be isotropic in the wind-wind collision zone. 

4.1. Inverse Compton scattering 

Inverse Compton (IC) scattering in the dense UV stellar radiation field of the massive 
main sequence star often dominates photon production rate at high energies. The compu- 
tation of the photon emission is based on the full Klein-Nishina cross section while for the 
IC losses of the electrons either the cross section in the Thomson limit or the Klein-Nishina 
approximations Eq. (9, 11, 14) are applied assuming the losses to be continuous. We have 
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further neglected triplet pair production, since the value of €tE /m e c 2 that we consider does 
not exceed 10 3 (Mastichiadis 1991). We restrict this work to long-period binary systems, for 
which the wind momentum from the WR-star clearly dominates, thus placing the WR-star 
at a large distance from the collision region. This allows us to neglect the stellar radiation 
field of the WR-star as a target photon field for IC scattering. For the stellar radiation field 
of the OB-star the monochromatic approximation 



with ex = 2.7/ceTgfj is employed. All seed photons are approaching the emission region 
from the same direction. In this case the full angular dependence of the IC scattering rate 
has to be taken into account, since the scattered power per volume element depends on the 
scattering angle # p h- In appendix C we calculate the IC photon production rate, n(e s , Q s ), for 
an arbitrary target photon field n(e, /i p h) that scatters off an isotropic electron distribution. 
We find a declining scattering rate with decreasing scattering angle 6* p h, in agreement with 
earlier works (Reynolds 1982a; Brunetti 2000, ; see Fig. C.l). The volume-integrated emitted 
photon power is calculated by 



where the integrals have been solved numerically. The IC flux variations directly translate 
into a change of IC power and maximum energy with viewing angle # p h (see Fig. 10). For 
a given system inclination the total emitted power and maximal scattered energy therefore 
varies with orbital phase (Fig. 11). These anisotropy effects may be detectable with near 
future 7-ray experiments like GLAST (see Sect. 6). For a non-negligible size of the convection 
zone the volume integration may lead to photon spectra that show a kink. This feature occurs 
at energies above which the convection zone is lacking high energy particles. The combined 
effect of both, a deficit of high energy particles in the convection zone and a visible increase 
of the total flux from the convection zone, results in the kink at 1-10 MeV in Figs. 10,11. 



Non-thermal relativistic bremsstrahlung losses are non-negligible in the stagnation point 
of the colliding winds where the (compressed) gas density may reach values of 10 5 ' 9 cm -3 
whereas in the convection region the decreasing gas density makes its contribution minor. 
The nonthermal bremsstrahlung photon flux using a typical ISM metallicity (90% H, 10% 
He, neglecting contributions from higher atomic number particles) 



n(e, /i ph ) = n S(e - e T ) - fx ph ) 



(28) 




(29) 



4.2. Relativistic bremsstrahlung 




(30) 
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is calculated numerically in the relativistic limit (valid for E/m e c 2 > 15/ Z, Z is the atomic 
number) using the differential cross section from (Blumenthal & Gould 1970). Examples 
of bremsstrahlung spectra are shown in Fig. 12 for the electron spectra in Fig. 3. The 
shape of the bremsstrahlung photon spectrum in the relativistic regime reflects the shape 
of the electron spectrum. The larger the binary separation is, the larger is the difference 
between the turnover energy from the Thomson to KN-loss regime and the cutoff energy 
(see Fig. 3). This leads to a decline of the bremsstrahlung spectrum at an energy that 
increases with the binary separation. An estimate of the thickness d of the emission region 
was provided by (Eichler & Usov 1993) who showed that d « x — x b- With r w 10 12 cm 
the total emitting volume in the acceleration region is nr^d 10 38 - 39 cm 3 . Fig. 12 shows the 
resulting relativistic bremsstrahlung power spectra for various binary separations. Despite a 
somewhat larger emitting volume for colliding wind systems with a large binary separation, 
the total bremsstrahlung emission declines with binary separation due to a rapidly decreasing 
target gas density in the collision region. For D = 10 14 cm the compound acceleration 
and convection region spectrum is also presented in Fig. 12. Due to the volume effect the 
dominant contribution to the total emission spectrum comes from the convection region at 
large r, thus increasing the overall bremsstrahlung intensity. At the same time the deficit of 
high energy particles in the convection region at large distances from the acceleration region 
(see Fig. 6) leads to a deficit of high energy photons. This causes the feature at ~ 100 MeV 
in the total bremsstrahlung spectrum shown in Fig. 12 (plotted for a D = 10 14 cm binary 
separation). 



4.3. 7r°-decay 7-ray emission 

Collisions between cosmic ray protons and nucleons in massive star winds are rather 
rare and occur on average a few times per year for gas densities in the wind collision region 
typical for long-period binaries like WR 140. 

The stationary proton spectra as shown in Fig. (8, 9) are used to calculate the 7r°-decay 
photon spectrum. Since the proton spectra above the threshold for pp-interactions in general 
reflects the shape of the acceleration spectrum, one expects pure power law particle spectra 
in this energy range. The formalism developed by (Pfrommer & Enfilin 2004) for the 7r°- 
decay 7-ray production of pure power law particle spectra seems therefore appropriate to 
use. The resulting 7r°-decay 7-ray spectra (calculated for a 4 He mass fraction of 0.3 for the 
wind metallicity; see Sect. 4.2) from the acceleration region are shown in Fig. 13 for various 
binary separations. The uppermost curve corresponds to the combined acceleration and 
convection region 7r°-decay spectrum for a binary separation of D — 10 14 cm. 



For typical parameters of colliding massive wind systems and maximal allowed injection 
power the radiative luminosity from 7r°-decay lies therefore typically on a < 10 33 erg/s flux 
level, leading to a 7r°-decay luminosity from wide binary systems that can in general be 
neglected when compared to the expected IC-luminosity >1 GeV, provided the emitting 
electron spectrum extends to > 10 4 MeV. 



4.4. 7 — 7 opacity in the stellar radiation fields 

Above ~ 100 GeV the optical depth r 77 due to photon-photon pair production (77 — > 
e + e~) in the intense stellar radiation field of the main sequence star may reach non-negligible 
values depending on its spectral type. This may modify the 7-ray spectrum escaping from 
the source by a factor exp(— r 77 ). 

The 7 — 7 opacity dr 77 along a path element dy in a radiation field with differential 
photon number density n(e, Q) is given by (Gould & Schreder 1967) 

dr 11 (E 1 , e, Q) = dyn(e, f2) <7 77 (e, E y , /i 77 ) (1 — /i 77 ). (31) 

With the target photon density 

r 2 

n(e) = n — s - — 5(e - e T ) 5{n - /x 77 ) 

x + y 2 + 2yfi ph Vx 

above the stellar radius r s , and x = x 2 + r 2 , this can be re- written as 

(t? \ ^o(e = e T )r2 f f°° sa 77 (s) 

r 77 (£ 7 ,/i ph ,;r,r) = / # / dy — — — H{s 

87r^ 7 e T y J x + y 2 + 2 Vxy/j, ph (9 L , <f) L ,cos9 = x/Vx, 0) 

(32) 

for photons of energy E 1 escaping from the emission region at location r along a path y with 
angle # p h = cos#l cos# + sin#£ sin# cos(0 — 0^). Here, s = 2_E 7 e(l — /i 77 (y)) is the square 
of the centre of momentum (CM) energy, cr 77 (s) the cross section for photon-photon pair 
production (Gould & Schreder 1967), s tiiI = (2m e c 2 ) 2 is the squared minimum CM energy 
required to overcome the process threshold, /i 77 (y) = cos^ 77 (|/) is the cosine of the angle 
between the two interacting photons, and cos# = x/y/x 2 + r 2 . 

In Fig. (14-16) we show the absorption optical depth due to photon-photon collisions in 
a T = 43000 K radiation field (i.e. e T = 10 eV) for different r OB = 1-2 - 24 • 10 13 cm being 
the separation of the emission region from the stellar photosphere, and for various angles 
9l = 0° — 180°, respectively. Large viewing angles decrease the process' threshold energy 
in the observer frame, and increase the opacity at the same time due to a longer path way 
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and the rise of the CM energy. Fig. 14 shows that variations of the optical depth with angle 
L by several orders of magnitude are possible, and absorption can be quite severe (up to 
t 77 ~ 100 for D = 10 14 cm). Except at large 8l the effect of varying r mainly impacts the 
absorption depth near threshold (see Fig. 15). The opacity is also strongly dependent on 
the radiation field density which is dependent on the massive star's luminosity as well as on 
tob- This can be observed in Fig. 16 where r 77 ranges between 0.002 to 0.04 at the peak of 
the cross section and along the 9 L — 45° sight line. 

This essay shows that 7-ray absorption due to pair production can not be neglected in 
general for massive star systems, but must be treated individually for each system and is 
strongly dependent on the sight line. 

5. Applications 
5.1. WR 140 

WR 140 is the archetype of a spectroscopic long-period massive binary system of col- 
liding winds which shows periodic dust formation around periastron passage as well as 
nonthermal radio emission at phases where the colliding wind region is outside the pho- 
tosphere for free-free absorption. The system has been well monitored at radio, IR, op- 
tical, UV and X-rays (e.g. Williams et al. 1990; White & Becker 1995; Panov et al. 2000; 
Setia Gunawan et al. 2001b). It consists of a WC 7pd and a 04-5 V companian (of bolo- 
metric luminosity Lboi = 1O 5 8 L and effective temperature T e ff = 43000 K) with an orbital 
period of 2899 ± 10 days. Located at a distance of 1.85 kpc (Dougherty et al. 2005) the 
large inclination system (i = 122° ± 5°; Dougherty et al. 2005) possesses an excentricity 
of 0.881 ± 0.04. With periastron passage being defined by phase $ = the argument 
of periastron is 00 = 47° (Marchenko et al. 2003). In this work we use the basic parame- 
ters that has been recently re-determined by (Dougherty et al. 2005). The mass loss rate of 
8.7- 1O~ 6 M / yr and the wind velocity of the O-star t>oB = 3100 km/s as compared to the mass 
loss rate of 4.3 • 10~ 5 M o /yr of the WR-star at velocity i>wr = 2860 km/s places the collision 
region at a distance of x ~ 0.32.D ps 9.6-10 12 . . . 1.5-10 14 cm from the O-star due to momentum 
balance and taking also into account the excentricity of the system. This is a factor ~ 2.5 
smaller than used in earlier works (e.g. Eichler & Usov 1993; Benaglia & Romero 2003). At 
this distance the wind gas density lies between ~ 4 • 10 6 . . . 9 • 10 8 cm~ 3 . For the strength of 
the shock we assume a compression factor c r = 4. The radio emission reaches its maximum 
roughly at phase 0.83. The nonthermal radio spectrum at phases $ = 0.671 . . . 0.955 shows 
typically a a = —0.5 spectrum, indicating that the underlying particle distribution obeys 
a power law with spectral index p = 2 (Williams et al. 1990; Dougherty et al. 2005). The 
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spectral requirements are in agreement with a convection velocity of ~ 1410 . . . 1540 km/s 
in the framework of our simplified geometrical picture. For predicting WR 140's high en- 
ergy emission from primary electrons and protons we have selected four orbital phases ap- 
proximately equally separated where its non-thermal radio emission has been determined 
(Dougherty et al. 2005), namely $ = 0.995, 0.2, 0.671 and 0.8. In order to reproduce the 
flux level of the cm radio data for those phases we require a energy injection rate of thermal 
particles that varies with orbital phase and lies between io 32 - 5 - 33 - 5 erg/s or 10~ 3 . . . 10~ 2 % 
of the OB-wind kinetic energy. The total emission volume V = r^^nd with d ~ x is 
also assumed to vary with orbital phase (we used r max = 4 • 10 14 cm for $ = 0.2, 0.671 
and 0.8, and r max = 10 13 cm for $ = 0.995) as suggested by the cm radio observations 
(Dougherty et al. 2005). The magnetic field at the O-star surface is assumed of order 100 G, 
which leads to field values of 0.2 . . . 3.5 G in the wind collision zone (in agreement with 
estimates for the equipartition magnetic field: Benaglia & Romero 2003). The diffusion co- 
efficient n a determines the acceleration rate and, together with the energy loss channels, the 
maximum particle energy. In order to allow for relativistic electron energies, necessary to 
explain the observed synchrotron radiation, k & must be low enough that acceleration gains 
are able to overcome the Coulomb losses. On the other side, Bohm diffusion is the limit for 
the diffusion coefficient k<±, which is connected to K a . This leads to a rather narrow range 
of possible values for the diffusion coefficients if one requires the production of > 100 MeV 
photons via the IC process at least at phases close to apastron. For our models presented in 
Fig. 18, 19, 20 we used K a = 3.76 • 10 19 cm 2 /s which results in an acceleration region size of 
r = 2.4 . . . 2.6 • 10 11 cm and escape time of T = 1.6 . . . 1.9ks. Detections of the maximum 
photon energy will shed light on the exact value of n a . 

Fig. 17 shows the steady-state electron spectra at orbital phases 0.955 (close to pe- 
riastron), 0.2, 0.671 (close to apastron) and 0.8 for these parameter values. Due to the 
excentricity of the system the wind collision region is located closest to the main sequence 
star at periastron while it is furthest away at apastron. All spectra except at periastron ex- 
tend to very high energies with the cutoff due to the finite size of the acceleration region, and 
a smooth change in the spectral index where Klein-Nishina effects set in. The sharpness of 
the cutoff is a result of the approximations employed in this work, and may appear somewhat 
softened in more sophisticated (however then indispensably non-analytical) calculations. 

At periastron the wind collision region is closest to the OB-star, and the electron spec- 
trum is cutoff already at ~100 MeV due to (Thomson regime) IC losses in the intense 
photospheric radiation field. Fig. 19 and 20 shows the corresponding broad band SEDs in 
comparison with the sensitivity of the INTEGRAL instruments, GLAST, AGILE and the 
current generation IACTs like H.E.S.S., MAGIC and VERITAS. The EGRET upper limit 
has been derived from the summed P1-P4 EGRET data at the location of WR 140 (see also 
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Miicke & Pohl 2002). At all orbital phases IC emission dominates over all other radiation 
processes except close to periastron. Here the IC spectrum cuts off already at a few MeV, 
corresponding to the cutoff in the electron spectrum. The 7-ray domain is therefore only 
covered by relativistic bremsstrahlung radiation (up to ~100 MeV following the electron 
spectrum) and 7-rays from the 7r°-decay. (Ion-electron) bremsstrahlung emission turns out 
to lie always below the IC photon output at a level that is not detectable for current and 
near future instruments. Assuming the maximal possible injection rate of thermal particles 
into the acceleration region, 7r°-decay 7-rays may possibly be detected with GLAST at or- 
bital phases close to periastron where the density of the target material is enhanced. At 
apastron the losses in the photospheric radiation field as well as due to the wind particle 
density are low enough to allow the photon spectra to extend into the 10-100 GeV regime 
with a IC flux level detectable, particularly with low-energy threshold IACTs (e.g. MAGIC) 
and GLAST. This is true despite photon absorption from pair production which leaves its 
fingerprint above ~100 GeV (see Fig. 20: at $ = 0.67 r 77 (r = r ,E lOOGeV) ~ 1 while 
at $ = 0.95 r 77 (r = r , E ^ lOOGeV) ~ 3). 

Fig. 18 zooms onto the IC spectra. The flux variations due to the anisotropy effect at 
different orbital phases are blurred by flux variations due to the changes in the radiation field 
density in the strongly excentric system. The feature at < 1 MeV stems from the deficit 
of high energy particles in the convection zone. For this parameter setting we predict a 
maximum of the IC flux at phase $ ~ 0.95, the minimum IC flux level would occur at phase 
~ 0.01. Thus the next maximum 7-ray flux level is expected around August 2008, just in 
time for a clear detection with GLAST. 



5.2. WR 147 

WR 147, a WN8(h) plus B0.5 V (with bolometric luminosity of L hol = 5 • 1O 4 L and 
effective temperature T eff = 28500 K, i.e. e T ~ 6.6 eV) massive binary system is among 
the closest and brightest systems that show non-thermal radio emission in the cm band. 
Owing to its proximity this system has been resolved in a northern non-thermal component 
(WR 147N) and a southern thermal one (WR 147S) with a separation of 575 ± 15 mas using 
the MERLIN instrument (Churchwell et al. 1992; Williams et al. 1997). The observed radio 
morphology and spectrum supports a colliding wind scenario for WR 147, as first proposed by 
(Williams et al. 1997). At a distance of 650 pc the implied binary separation is estimated to 
417 AU. The mass loss rates (M W r = 2.5-lO~ 5 M /yr, M OB = 4-lO~ 7 M /yr) and wind veloc- 
ities (%r = 950 km/s, vob = 800 km/s) place the stagnation point at 6.6 • 10 14 cm, in agree- 
ment with the MERLIN observations. A comprehensive study of the WR 147's radio emission 
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and the geometry of the system has been presented by (Setia Gunawan et al. 2001a). Any 
excentricity nor the inclination of the system are known so far, hence we assume % = 90° and 
e = for this application to WR 147. The non-thermal flux component can be well fitted by 
a power law with spectral index a = —0.43 with, however, poor statistical significance. For 
our presentation we use again the canonical value of a — —0.5 corresponding to an E~ 2 power 
law particle spectrum. Assuming strong shocks, a convection velocity of 399 km/s and a 
B-star surface magnetic field of 30 G (translating into a 25 mG field strength in the wind col- 
lision region, close to its equipartition value (Benaglia & Romero 2003)) the observed radio 
flux level can be reproduced if ~ 0.15% of the OB- wind kinetic energy is injected as thermal 
particles into the emission region of a total estimated volume of ~ 2.6 • 10 5 AU 3 . Again the 
diffusion coefficient is chosen such to overcome Coulomb losses and allow the production of 
> 100 MeV through the inverse Compton process without violating the Bohm-limit at the 
same time. With K a — 2- 10 21 cm 2 /s convection is dominant in ~ 90% of the total emission 
region, and the acceleration site covers a size of ro = 5 • 10 13 cm (escape time is To ~ 10 6 s). 

In Fig. 21 the resulting steady-state electron spectrum is shown. For a negligible system 
eccentricity it is the same at all orbital phases. The steepening at 10 4 "' 5 MeV is due to 
synchrotron losses, the slight upturn towards lower energies reflects the influence of the 
Coulomb-losses. Above ~ 10 5 5 MeV the Bohm-limit cuts off the particle spectrum. Fig. 22 
demonstrates the orbital IC flux variations due to the anisotropic nature of the IC scattering 
of more than one order of magnitude for the chosen system inclination. The maximum flux 
level is expected when the WR-star is behind the OB-star along the sight line. These flux 
changes will be in principle detectable for GLAST at all orbital phases, and also for low- 
threshold IACTs provided high particle energies are reached, which by itself depends on the 
acceleration rate. Thus measurements with IACTs will be able to place important constraints 
on the acceleration efficiency in these environments. The EGRET upper limit is taken from 
(Benaglia & Romero 2003). For an estimated wind gas density of up to 1.4 • 10 5 cm -3 non- 
thermal bremsstrahlung emission stays several orders of magnitude below the IC flux level, 
not detectable with current instruments to date. Also the contribution from 7r°-decay 7-rays 
remains negligible, even for a maximal possible thermal particle injection rate. For WR 147, 
being a extremely long-period binary with an assumed small excentricity, the photospheric 
UV radiation field at the location of the collision region is low enough to allow here the 
neglection of 7-ray absorption due to photon-photon collisions. 
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6. Detectability from soft to very high energy gamma-rays 

6.1. Soft gamma-ray instruments 

The regime of soft gamma-rays, here considered to be typical ranging from 511 keV 
up to several MeV, is currently only accessible by the IBIS and SPI instruments aboard 
INTEGRAL. Previous mission like SIGMA/GRANAT and COMPTEL did not reported 
the detection of soft gamma-ray emission from the WR systems studied here. Given the 
continuum sensitivity of the IBIS and SPI instruments 1 , accordingly scaled for a common 
observation time of 10 6 sec in order to achieve a 3 a detection, the detection of WR 140 
appears to be possible at soft gamma-rays only if significantly more observation time than 1 
Msec will be dedicated to observations towards such an object. Hard X-ray emission may be 
seen by ISGRI, most favorably in orbital phases where the line of sight towards an observer 
is parallel to the contact discontinuity of the wind collision zone or when both the WR and 
O star are nearly aligned and the wind collision zone is most pronouncedly exposed towards 
the observer. At soft gamma-rays, even in the most favorable orbital states, several Msec 
may be required to pick up emission from WR 140 (Fig. 18). This situation is very similar 
for WR147 (Fig. 22). 

6.2. High-energy gamma-ray instruments 

At energies between 30 MeV and 10 GeV, detection claims of gamma-ray emission from 
WR binary systems have been made already from COS-B observations. More particular, 
WR 140 (HD193793) was considered to be associated with the variable COS-B source 083+03 
(Pollock 1987), a source exhibiting a photon flux at the 5-10~ 7 ph cm _2 s _1 level at E > 300 
MeV. However, this individual association could not be confirmed by observations of the 
EGRET instrument aboard CGRO, which would have seen a source at the ~ 1.5-10~ 6 ph 
cm _2 s _1 at E > 100 MeV. Although a number of positional coincidences between unidentified 
EGRET sources and colliding wind binary systems has been noted by (Kaul & Mitra 1997) 
and (Romero et al. 1999), the individual case for an association between WR 140 and the 
unidentified EGRET source 3EG J2022+4317 appears to be vague on the basis of the given 
observational evidence at 7-rays. 3EG J2022+4317 is cataloged (Hartman et al. 1999) with 
a flux of 2.5-10~ 7 ph cm~ 2 s~ 1 , and the spectrum fitted with a power-law with an index of 2.3 
± 0.2; it is further to be characterized by a rather irregular uncertainty contour pointing 



INTEGRAL AO-3 documentation, SPI and IBIS Observer's Manual, ESA 2004, 
http: / / www.rssd.esa.int /Integral/ A03/ 
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towards source extension, and indication of source confusion above E > 100 MeV. The 
EGRET observations were taken within the first four years of the CGRO mission between 
1991 and 1994, with the gross obtained toward the periastron phase of the binary system. In 
the sparse EGRET observations above the detection threshold of the 3EG catalog there is 
no evidence for variable 7-ray emission (Nolan et al. 2003). Since WR 140 is at a distance of 
0.67° to the nominal position of 3EG J2022+4317, and only barely consistent with the 99% 
source location uncertainty contour, which prevents any conclusive identification between 
both objects, it is appropriate to determine an upper limit at the position of WR 140 under 
consideration that 7-ray emission of 3EG J2022+4317 is fully taken into account but not 
been associated with WR 140 (Miicke & Pohl 2002). This upper limit is consistent with our 
predictions given in Fig. 18, considering that EGRET observations were performed over a 
superposition of orbital states including the periastron phase, where the cutoff due to IC 
losses in the Thomson regime will not allow any detectable 7-ray emission at all at energies 
above a few MeV. Since AGILE will exhibit a similar sensitivity characteristics compared to 
the EGRET instrument, chances for AGILE to clarify on 3EG J2022+4317 are only given 
if AGILE is operated in orbit over a long period. The instrumental sensitivity of AGILE 2 
is indicated for a 5 a detection on the basis of 10 6 sec observation time. Thus, it'll be the 
Gamma-ray Large Area Space Telescope (GLAST) to give an observational reassessment of 
the unidentified EGRET source 3EG J2022+4317 and its association or non-association with 
WR 140. Clearly, GLAST 3 exhibits a sensitivity characteristics to achieve a detection not 
only from observations accumulated over various orbital states of the colliding wind binary, 
it has a dedicated chance to provide results from individually selected orbital states. As 
mentioned above, during periastron phase there is no high-energy 7-ray emission predicted 
(Fig. 19), but the changes in the 7-ray flux of WR 140 when going into periastron or coming 
out of periastron phase will be detectable, since the sensitivity of GLAST will be indeed 
sufficient to detect this system in the more favorable orbital states when the edges of the 
contact discontinuity of the wind collision zone lines up in the line-of-sight or the wind 
collision zone is most extensively exposed towards the observer. GLAST will have this 
detection potential up to energies where the the finite acceleration site cuts off the emission, 
approximately up to ~ 50 GeV. The chance to detect high-energy emission in case of WR 147 
with future instruments is even somewhat better than for WR 140. If AGILE with its better 
instrumental point spread function than EGRET will be able to distinguish the location 
of WR 147 from the bright 7-ray source 3EG J2033+4118, its sensitivity will be sufficient 
to observe 7-ray emission over the majority of favorable orbital states (Fig. 22). However, 



2 lit t p : / / agile . mi . iasf . enr . it /Homepage / performances . shtml 
3 http://www-glast. slac.stanford.edu/software/IS/glast Jat_performance.htm 
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the periastron phase will still be out of reach for AGILE, see Fig. 23. GLAST, again, 
will have the instrumental capability to both distinguish a test position at the location of 
WR 147 (I = 79.85,6 = -0.32) from the presumably bright 7-ray source 3EG J2033+4118. 
Depending on a possible contribution of further sub-threshold sources in the vicinity of 
3EG J20233+4118, this might not be too straight forward, though. 

6.3. VHE gamma-ray instruments 

Complementary to the satellite experiments operated at soft and high-energy 7-rays, 
the continuous development of the Imaging Atmospheric Cherenkov Telescopes (IACTs) to- 
wards better sensitivity, even more improved angular resolution, and lower energetic thresh- 
olds, clearly exhibit the capability to detect 7-ray emission produced in colliding winds of 
massive stars in binary systems. Especially the chance to accumulate a wealth of photons 
on sub-hour time scales will enable the IACTs to test the predicted orbital flux variations 
precisely. A unique chance to detect emission from WR 140 is given for those experiments 
able to work at the lowest possible threshold, i.e. MAGIC below 100 GeV. Depending on 
the actual shape of the cutoff due to the finite site of the acceleration site, events may be 
seen only in extreme low-energetic event selections, and the respective sky location is fur- 
ther characterized by an absence of higher energetic photons. For an array of IACTs like 
VERITAS (Weekes et al. 2002) the best sensitivity is achieved at energies above the cutoff 
due to the finite acceleration site in WR 140, subsequently making any detection prospects 
heavily dependent on the actual shape of the cutoff. IACT arrays located in the southern 
hemisphere like H.E.S.S., or CANGAROO may not have the chance detect this system at 
all due to the higher energetic threshold when observing under low zenith angle conditions. 
In case of WR 147, where the orbital parameters are more promising to detect high-energy 
7-ray emission up to several hundreds of GeVs, low threshold IACTs located in the northern 
hemisphere will have a distinct chance to detect WR147 in non-periastron orbital phases 
(Fig. 22), provided that the 7-ray analysis will not introduce hard cuts for preference of a 
higher energetic event selection. For both satellite and ground-based instruments the collid- 
ing wind zone will not appear to be spatially resolved, presenting individual colliding wind 
binary systems as point-source candidates at the 7-ray sky. 

7. Conclusions and discussion 

In this work we have calculated the emission from non-thermal steady-state particle 
spectra built up in the regions of colliding hypersonic winds (assumed to be homogeneous) 
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of long-period massive binary systems with the stagnation point defined by balancing the 
wind momenta and under the assumption of spherical winds. The shocked high-speed winds 
are creating a region of hot gas that is separated by a contact discontinuity. The gas flow 
in this region away from the stagnation point will be some fraction of the wind velocity 
which we kept constant here. A simplification of the geometry from a bow-shaped to a 
cylinder-shaped collision region allowed us to solve the relevant diffusion-loss equations ana- 
lytically. We considered first order Fermi acceleration out of a pool of thermal particles, and 
took into account radiative losses (synchrotron, inverse Compton including Klein-Nishina 
effects, bremsstrahlung and Coulomb losses), (energy- independent) diffusion by introducing 
a constant escape time T and convection/advection with constant speed. Above a certain 
distance from the stagnation point convection dominates over diffusion with the transition 
point determined by balancing the diffusion and convection loss time. Correspondingly, we 
devided the emission region into a region where acceleration/diffusion dominates, the "accel- 
eration zone", and the outer region where convection/advection dominates, the "convection 
zone" . 

Electrons may reach relativistic energies, once they overcome the heavy Coulomb losses 
in the dense shocked material, through diffusive shock acceleration up to the Bohm diffusion 
limit. For wide binary systems this latter constraint is often severe, while for close binaries, 
radiative losses mostly cause the cutoff. Taking into account existing upper limits of the 
stellar surface magnetic field strength of massive stars, inverse Compton losses in general 
dominate over synchrotron losses if in the Thomson loss regime. We have shown, however, 
that losses may well extend into the transition region leading to the extreme Klein-Nishina 
regime. The flattening of the Compton loss rate there may in some cases cause the syn- 
chrotron losses dominate eventually. Thus a rigorous treatment of the Compton losses must 
include Klein-Nishina effects. For this purpose we have derived analytical approximations 
for Klein-Nishina losses that are suitable for massive binary systems, and allow at the same 
time to solve the relevant diffusion-loss equation analytically. Despite the high density en- 
vironment of the emitting collision region, non-thermal bremsstrahlung losses prove general 
to be of minor importance. This turns out to be true also for the corresponding radiation. 

We have studied inverse Compton radiation, the main emission channel for relativistic 
electrons in these systems, in severe detail. The use of the full Klein-Nishina cross section 
leads to a spectral softening at the high-energy end of the emitted radiation. Since the stellar 
target photons for inverse Compton scattering arrive at the collision region from a prefered 
direction, the full angular dependence of the scattering process has to be considered. Its 
anisotropic nature leads to variations of the flux level by up to several orders of magnitude 
(depending on system inclination and eccentricity) as well as cutoff energy with orbital phase. 
The maximum flux and cutoff energy occurs when the WR-star lies behind the OB-star. We 



-27- 



consider therefore massive binary systems as 7-ray sources that are variable on the time scale 
of their orbital period even in the absence of a strong system eccentricity. The inclusion of 
convection/advection effects into the calculation of the particle spectra reveals a possibly 
visible spectral feature, too. Because of a deficit of low-energy particles in the convection 
zone, a softening of the volume-integrated radiation spectrum may occur if the convection 
zone is sufficiently large compared to the acceleration zone. A detection of this feature would 
give valuable information about the particle propagation properties in the emission region. 

Since thermal protons are most likely wind constituents as well, diffusive shock acceler- 
ation implies the presence of relativistic protons in the wind collision region. If they reach 
energies of several GeV, their presence may show up as 7r°-decay 7-rays produced through 
inelastic proton-proton collisions. Their detectability, however, depends not only on the rel- 
ativistic electron-to-proton ratio and the instrument capabilities, but also on the importance 
of the competing radiation mechanisms. E.g. in the case of WR 140 close to periastron, 
the otherwise dominating inverse Compton radiation most likely cease to reach sufficient 
high energies that would allow MeV-GeV emission, which increases the chance of detecting 
7r°-decay 7-rays. 

Finally, we find that photon-photon pair production can not be neglected if the produced 
radiation exceeds energies of ~ (kT/eV)' 1 TeV, which lies typically at 50-100 GeV. The 
absorption optical depth thereby depends sensitively on orbital phase and system inclination. 

Although many free parameters are involved in the presented model for high energy 
emission from the wind collision region of massive binaries, few are those which are unre- 
lated to observations, and even fewer those which - if changing - may have a significant 
impact on the predicted 7-ray intensity. Indeed, since IC emission seems the dominant emis- 
sion process at high energies in most cases, the high energy output can directly be deduced 
from the knowledge of the synchrotron emission. While the non-thermal radio flux level pro- 
vides information on the required injected energy in form of electrons if the magnetic field 
is known, its radio spectrum constrains acceleration and propagation parameters. Equipar- 
tition arguments together with lower magnetic field limits from the Razin-Tsytovich effect 
(e.g. Chen 1992; Benaglia & Romero 2003), and observational limits on the stellar surface 
magnetic field in connection with a plausible dipole field configuration can be used to esti- 
mate the magnetic field strength in the collision region. The fact that relativistic electrons 
exist, supplies a lower limit on the acceleration rate, while an upper limit is given by the 
Bohm diffusion regime for the likely case of diffusive shock acceleration operating in these 
objects. 

We applied our model to two archetypical WR-systems: WR 140 is arguable the most 
popular among these sources. We predict WR 140 to be detectable with GLAST and 
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MAGIC, except at phases close to periastron due to an early cutoff of the electron spetrum 
already at ~ 100 MeV. This may lead to the dominance of bremsstrahlung and hadronically 
produced 7-rays above ~ lMeV, at this phase, while at phases far from periastron inverse 
Compton radiation is predicted to dominate at all energies. Orbital flux variations at high 
energies far from periastron are expected with amplitudes that vary by a factor ~ 2. 

The ~factor 10 wider binary system WR 147 is notable for being the brightest (because 
closest) system at radio frequencies, and for being one of the few systems where the thermal 
and non-thermal radio emission are observed to arise from spatially different resolved regions. 
Due to a lack of knowledge of the system parameters we model WR 147 face-on with no 
significant eccentricity. The low target photon density at the collision location makes photon 
absorption negligible here, and at the same time allows the electron spectrum to extend up to 
sub-TeV-energies if the acceleration efficiency is favorable. This would lead to radiation up 
to the 100 GeV-region on a flux level possibly detectable even with VERITAS at some orbital 
phases, while GLAST has good chances to trace this system at all phases. INTEGRAL'S 
sensitivity at 7-ray energies will most likely be insufficient to discover these sources as 7-ray 
emitter given the finite amount of observation time in individual instrumental pointings. 

In this work we concentrated on the emission from a steady-state particle spectrum. 
Being time-dependent systems in general, a time- dependent diffusion-loss equation shall give 
a more realistic description of the emitted intensity. In this case we expect that, similar to 
supernova remnants, the electron spectrum will slowly be built up with the maximum particle 
energy increasing with time. Typically the electron spectrum is fully developped after a few 
tens of hours. The uncertainty of a given phase corresponding to the here calculated steady- 
state emission therefore lies typically in this time range. A more comprehensive discussion 
of the behaviour of massive colliding wind systems in the framework of a time-dependent 
diffusion-loss equation will be considered in a forthcoming paper. 

In conclusion, we consider colliding wind regions of massive binary systems that are wide 
enough to avoid radiative braking, as promising sources of high energy emission that may 
extend far beyond the X-ray band. High energy observations of these systems by sensitive, 
low-threshold IACTs and satellite instruments can be used not only to derive geometrical 
details but also to explore the efficiency of diffusive shock acceleration at densities much 
higher than in other astronomical objects with high Mach number shocks, e.g. supernovae. 
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A. Solutions for steady-state electron distributions in the acceleration region 



Here we derive the analytical solutions of Eq. 4 for electrons suffering radiative losses 
following Eq. 8. 
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It is an exact solution for E < E TL . With the approximation Eq. 9,11 for the inverse 
Compton losses Eq. (A2, A4) describes the electron spectrum for E > E s if one substitutes 
^syn&ic — ¥ b sy n, ^br - ► b hr + q a and 6 cou i — > 6 cou i + <2V In practice, however, relativistic 
bremsstrahlung and Coulomb losses above E > -&tl can usually be neglected in the systems 
considered here. For _&tl < E < E s the solution of Eq. 4 reads: 
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This preprint was prepared with the AAS IATgX macros v5.2. 



-33- 



( E%a - b hr ) - b synicIC E + b lc , TL Ej/E g ) \ -VP<-W*b] 

V E 2 ((a - b hl ) - b syDgzlc E + b lc ^E*/E g ) ) 1 ) 



exp 



frsyn&IC / 26 ICjTL /E 9V /|A|(E - E 



(a-6 br )T 0V ^A] V (26i C ,TL^/^-6syn&ic)(2&ic,TL^o/^-&s y n&ic) - A 



B. Solutions for steady-state nucleon distributions in the acceleration region 

Here we derive analytical solutions for steady-state nucleon spectra solving Eq. 4. 
At non-relativistic energies and below the Coulomb-barrier the solution is: 

On ( E \ ~ («- b bel)T0 

for E < E c and N(E > E c ) = 0, while above the Coulomb-barrier E > E m we find: 

for E > E m , E c = (6 ab /a) 2 / 3 , 
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At relativistic energies > m p c 2 = E re \ and above E > E t \ 1T , the pp-interaction threshold, 
one finds: 
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for E < E m , 
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for E le i < Eq < -Ethr- 

Finally, at relativistic energies but below E < E thl the steady-state spectrum turns out 
to follow: 
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for i?o < -E m , and 



for E m < E < E Te \. 



C. Analytical calculation of the inverse Compton scattering rate 

In the lab frame we choose a polar coordinate system such that the line-of-sight marks 
the z-axis. A single incident electron is then fully described by its Lorentz factor, 7, and the 
polar angle, \i e = cos# e , and azimuthal angle, <f> e , to mark its direction of flight. 

The photon field is described by the differential photon spectrum 

n ^ nph) = W^ h (C1) 

where e = E ph /m e c 2 is the dimensionless photon energy. 

The scattering rate (quantities of the scattered photon are indexed with s) for a given 
differential number density of electrons, n e (7,f2 e ), is 
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where the part in brackets is the scattering rate for a single electron and fj, p h, e is the cosine 
of the angle between the photon and electron flight directions, i.e. 



AW = VeVph + V 1 - A*l f4h COS & (C3) 

The differential cross section is well known in the electron rest frame. It is therefore attractive 
to calculate the scattering rate per single electron in its rest frame (ERF, indicated by an 
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asterisk) and then to transform the result back into the lab frame. Since then only photon 
spectra have to be transformed we can use the invariants 

7* 6 (C4) 



C.l. The scattering rate for a single electron 



In the following we will assume that the electrons have an isotropic distribution, in which 
case the scattering rate can not depend on the azimuthal angle of the incoming photons. 
Therefore we can set <p p h = 0. 



The Lorentz transformation relating lab frame quantities to those in the ERF are 
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We also need to know the photons azimuthal angle with respect to the electron. Let the 
azimuthal angle of the line-of-sight be 4>i os ,e = 0. Then 
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Note that 
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The differential cross section in the ERF is 
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where £ is the scattering angle. 
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Given the photon angles after scattering, fx* and (/>* = Sj6 , the scattering angle can 
be calculated as 

COs£ = /^X,e 



+ V 1 - AW V 1 - ^/ cos(0 pV - s , e ) (CIO) 

where 

COs((t) p h,e - <t>s,e) = COS <fiph,e COS (fis,e + Sm (f) ph , e Sin SjC (Cll) 

Then the scattering rate in the ERF is 

n*(e:,{i:)=cj f^- e n ph (e,Q ph )^^dQde (C12) 
which can be transformed back into the lab system using Eq. C4 and 

= ttm: (C13) 

£ «=« i+ *-)=^k) (C14) 

Finally we obtain the scattering rate in the lab frame for a population of electrons 

h(e s ,Q s ) = J <j) J (j> dedQphdjdQe 
c (1 - Ptip hte ) da 

7 (i - pil^) n ^ a » h) ne(7 ' fie) (C15) 



C.2. The scattering rate for an arbitrary isotropic electron spectrum 



Note that the line-of-sight is defined by (fi s e = and /i se — \i e . Therefore cos(<j) p h, 

7 2 (1-/8/0* 



S|C ) = cos0 p/lie . Also [1 - /i* 2 ] = -a -p , hence 



COS <^ = - 

(1 -Pfl e )(l -PHph,e) 

= i 1 ~ ^P h (C1F,) 

7 2 (l-^e)(l-^,e) 



1 - /ie/ipfe.e - y/1- fJ.lJl- fJ% h COS < 

1 - COS e = T7 r-r-t- r (C17) 

f l-fe l-^, e 
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e*(l — cos£) = e • 



1 - flet*ph,e ~ V 1 - t4h,e COS <i>ph, 



1 - Uph 

= e- 



7(l-/^e) >(l-^e) (C18) 

The scattering rate does not depend on <fr ph , this integral thus being trivial. For the angles 
e and fj, e one is trivially performed by using the delta-functional in the cross section and 
the other one has to be done explicitely. 



The delta functional in Eq. C9 can be rewritten as 

4 r - T (i-w +f (i-J = 5 

With n p h(fiph) = f d(p p h n p h(Q p h) the scattering rate is 

r^cn e (7, Q e ) n ph (e, /i ph ) (1 - (3fi ph ,e) 



(C19) 



de / c?/ip ft d'j (b dQ 



2[7(l-/fy e ) + e(l-AV0]' 



*[/(^e,A*pft>0e)] 



(C20) 



x 



1 - 



1 — /V» 



! (1-A*pfc) 5 



7 2 (1 - /9/i e ) (1 - (3fip h ,e) ) 7(1 - /3/i e ) [7(1 - P^e) + e(l - fi ph )\ 
Given the argument / of the delta-functional 



1 - P(J, p h,e 



and 



as well as 



7(1 - Pn e ) + e(l - Uph) 7e 
1 ? 



(1-Pfi phte ) 7(1 -/fyO + e(l ~/v0 



(C21) 



(C22) 



<*(/) = [7(1 - /3/i e ) + e (1 - /vO] ^ [g = e s (7(1 - /tye) + e (1 - n ph )) - £7(1 - /3// pft , e )] 

(C23) 

Using the isotropy of cosmic ray electrons, n e ( 7 , f2 e ) = -^n e {j) we then obtain for the 
scattering rate 



x 



e(l - // pfe ) 



7e 



1+1- 



+ 



%(A*e,AV»>&)] (C24) 



! (1 -^phf 



76 s (1 - /3fJL e ) [7(1 - /3^ e ) + e(l ~ /VO] / 7(1 - /fye) [7(1 ~ /tye) + e(l - /vO] 

Now we need to find the zeros of the delta functional. Inserting Eq. C3 into g = we obtain 



1 + -(1 - Uph) 

7 



-e - P(e s - /jLp h e) pL e = -f3ecos<f) e yjl - fi* h y/l - ji 2 e (C25) 
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which is of the form 

A-Bn e = -Cy/l^ (C26) 

where C can be positive and negative, depending on <p e . Also A > and B > for the 
interesting case e s > e. Real zeroes fx of the above equation exist, if B 2 + C 2 — A 2 > 0. 



C.3. The case /j, ph = ±1 

If n p h = ±1 we have C = and then /j, — ^. The delta-functional transforms as 

%) = ^(/<e-§) (C27) 

and 

A(£ „ n .) = ^[ d , [ dl b.(7)^.^ = ±1) e (B _ _ 4) (C28) 

e(l - fXph) \ 2 e 2 (l - fi ph ) 2 



x 



1+1- 



7e s (1 - /3 no) [7(1 - /3a*o) + e(l - /vJ] / 7(1 - /3/x ) [7(1 - /?A*o) + e(l - A*ph)] 



C.5J. ^h = +l 

Here obviously -B/A = /? < 1, which implies that the argument of the delta functional, g, 
has no zero in the range of integration, and thus the scattering rate in the forward direction 
is precisely zero. 



Here 



and then 



Then also 



C.3.2. fiph — —1 



B = f3(t s + t) A = e s (l + — )-e (C29) 

7 



He = I < 1 =► e s < ^Y-p I e ( C3 °) 

l-/3/i = -^ (C31) 
7 e s + e 
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and 



4e 2 



7(1 - M [ 7 (1 - M + e(l - li ph )] = [l 2 - (e s - e) 7 - ee s ] (C32) 



Therefore 



f de f dl n e ( 7 )r^(e,^ = -l) Q 



/3 7 2 e(e s + e) 



1 + 



7 M 



(C33) 



x 



1+1- 



e s + e 



2ee s [7 2 - (e s - 6)7 - ee s 
where 6 denotes a step function. 



+ 



[l 2 - (e a - 6)7 - ee s ] 



C.4. The case \i vh 7^ ±1 
We obtain formally by squaring Eq. C26 



AB ± \C\VB 2 + C 2 - A 2 
- b 2 + C 2 

Since we have squared the original equation, not all formal solutions may apply. The am- 
biguity in sign is combined with the modulus of C, therefore writing C instead of |C| does 
not introduce new solutions. Then, assuming A > and B > 0, i.e. upscattering of photons 
in energy, we find 

AB ± CVB 2 + C 2 -A 2 
»e = (C34) 

is a solution to the original equation 

case © if C < jVB 2 + C 2 - A 2 (C35) 

case e if C < -^Vfi 2 + C 2 - A 2 (C36) 
The delta functional in Eq. C24 transforms as 

B ^ ~ A s(fjL e - At ) (C37) 

y ' -fCVB 2 + C 2 -A 2 V ; V ; 

Then 

- / * /„ Ifi+fc*) (C38) 
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x ( 5 ^o - A) (±^B 2 + C 2 - A 2 - C^j 



x 



1+1 



:(1 - 



+ 



e 2 (l - // p/l ) 5 



7e s (1 - /3^ ) [7(1 - PfM>) + e (l - fhh)] J 7(1 - /5a*o) [7(1 - /?A*o) + e(l - /vO] 

The integration variable e appears only as the argument of a cosine function, which itself 
appears only in the term C. We can therefore substitute C for <f> e . With C\ = (3e^J 1 — [i 2 h 

Ci 



(C39) 



and we obtain 



r 2 ce s 
An 



de I dji. 



'•ph 



Cl d(J n e ( 7 ) np h (e, M (£ 2 + C 2 - A 2 ) 



x 



1+1- 



-d i 2 eCy/Cf - C 2 V^ 2 + C 2 - A 2 
x (5/i - A) (±^V£ 2 + C 2 - A 2 - C^j 

e 2 (l-^) 2 



(C40) 



e(l - //ph) 



7e s (1 - /3 no) [7(1 - Pfio) + e(l - /v>)] / 7(1 - (3fi ) [7(1 - l3fi ) + e(l - /vO] 
It is now useful to discriminate two cases. 

1-/3/V* _ 7 2 e(l+/3)(l-/?AV») 



B > A 



e,<e 



1- (3+^(1- fi ph ) 1 + 7 e(l + - Aiph) 
for which only the © solution exists for all C. The second case is 



with 



e s2 = 



B < A and C 2 > A 2 



67 [7(1 - + e(l - fi ph )} 

1 + 2e 7 (l -fi ph ) +e 2 (l -/iph) 2 



e sl < e s and e s < e. 



s2 



1- 



1 + 2e 7 (l - fjL ph ) + e 2 (l - /vjj 
[ 7 2 (1 -/?AV») + e 7(l -Vph)} 2 



(C41) 



(C42) 



(C43) 



for which both the © and the solutions exist, but only for negative C < —C 2 = 
— \J A 2 — B 2 . If fiph is not very close to 1, e s2 does in fact exceed e s2 , but only by a very small 
margin. The relative difference (e s2 — € sl )/e sl < 10 -2 for 7 = 10 and it decreases rapidly 
with increasing 7. 
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C.4-1- The case B > A 



The integration variable C is an argument of /xo- We may substitute /x = /xo for C and 

use 

- g*^ - 1 - 1 ( C44) 

J- Cl Cy/Cf - CWB 2 + C 2 -A 2 sfW+Clh- V 7 ^" - A*) 

with 



^ = ^ ■ ^ 1 (C45) 



Aff ± Ci ^JB 2 + C 2 - A 2 
B 2 + C? 

to obtain the scattering rate as 



4:71 



X 



or 



1+1- 



fde fd» ph fd>y Wl^'M Q{B _ A) f" + <fr 
J J J p'j 2 e v /e 2 -2fi ph ee s + e 2 ^(/x - /x_)(/x+ - /x) 

e(l - /x pft ) \ 2 e 2 (l - /x p/l ) 2 



7 e s (1 - (3fi) [ 7 (1 - /3/x) + e(l - /x pft )] 7 7 (1 - /?//) [7(1 - Pn) + e(l - /x p/l )] 



X 



47r y J P J P^e^/e 2 - 2fi ph ee s + e 2 J ^ ^(/x - /x_)(/x+ - /x) 

9 e 2 e s (l - xip fe ) 2 - 2e(l - /ipfc) e 2 (l-/x p/l ) 2 



7 e s (1 - /3/x) [7(1 - /3/x) + e(l - AvOl ' e 2 ( 7 (1 - /3/x) [ 7 (1 - /3/x) + e(l - av)]) 2 

(C46) 

which is essentially the sum of three integrals I\, I2, I3, which are defined by the terms in 
brackets. 

For the first integral note that 

I = = — arcsm I = n =>- l 1 = 2it (C47) 

h- vO" - A*-)(A*+ - A 1 ) V A*+ — A*- / 

The second integral has the structure 

/ 2 = f + ^ = 7 ^ (C48) 

where 

a=-I ^-i-^ (048) 
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We have 



h = D 



7/9 



e (! - fJ-ph) 



d/j, 



d\i 



Io——D 



1(3 



e(l - fiph) 



dt 



1 vM + 5 a t + C a t 2 7 1 + B b t + C 6 t 2 



(C50) 

^ dt 1 (C51) 



where with (a, 6) standing for a or 6 

-4a = -1 #a,& = //+ + //- + 2 (a, 6) 
C a ,b = - (/i- + H+){a, b) - (a, bf 

Because 

C a , 6 < and A = 4A a C a , b - 5 2 fe < 
the solution of the integrals in Eq. C51 is 



1 



dt 



IT 



M— +(a,b) 



hence 



h — — 7T 



ee s (l - fiph) - 2 



The third integral is of the form 



h = D 



a/ (/!+ + b) (JjZ + 6) a/(/U+ + a)(//_ + a) 



where 



We may write 



M _ y/(fj, - (J.-){fJ.+ -fJ.)(fJ- + a f(^ + bf 
e 2 (l- fi ph f 



D = 



1 



(/x + a) 2 (/i + 6) 2 (a-6) 5 



1 



1 2 



1 



so that 



(a - 6) 2 
/ 3 = 



1 1 

+ 



fi + b ji + a 
2 



(/i + 6) 2 Gu + a) 2 (/i + a)(/i + 6) 



h- y/(fJ. 



(C52) 
(C53) 

(C54) 
(C55) 



(C56) 



(C57) 



(C58) 



(C59) 
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x 



1 



1 



_(/i + 6) 2 (/i + a) 2 (fi + a)(fi + b) 
= Ji + J 2 + h 
The integral J 3 is formally identical to I 2 



J* 



2tt 



7/3e 2 e(l - fiph) 



1 



y(^+ + 6)(/i_ +6) v/(^+ + a)(/i_ +a) 
The other integrals are (for ,J\ replace a by b) 

n ji + + /!_ + 2(a, 6) 



^2^,2,02 t 

e s 7 P J 2 = - 



2[( / , + + (a,6))(^ + (a,6))] 3/2 
The total scattering rate is then 

h ( e s , Q 8 ) = T -^ ± J ' de J d/i ph J dry 
1 



nAj) n ph (e,n ph ) 



0(e sl - e s ) 



x 



/i + + yU_ + 26 /i + + /x_ + 2a 



4e2 7 2 /? 2 v>/^++6)^-+6) +«r 



(C60) 



(C61) 



(C62) 



(C63) 



27/3 lPe 2 s e(l - fi ph ) J \y/(fj, + + &)(//_ + 6) v /(// + + a)(/i_ + a)_ 

Note that — is an extremely small number as well as the difference of the. 

parameters a and b. For ji p h = — 1 we reproduce the earlier solution (C33). 
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Fig. 1. — Model geometry for the collision of two stellar winds of a WR+OB binary system. 
A strong shock is created in the wind collision region which is situated at a distance x = xob 
from the OB-star. This region of thickness d consists of an acceleration zone of size r$ where 
suprathermal particles from the stellar wind are accelerated, and an adjacent convection 
zone where particle streaming along the wind contact surface dominates over diffusion. For 
details see text. 
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Fig. 2. — Energy loss rates due to inverse Compton scattering (thin solid lines) in the 
Thomson regime (TL) and Klein-Nishina regime (KN), synchrotron radiation (dashed-triple- 
dotted line), relativistic electron-ion bremsstrahlung (dotted line) and Coulomb interactions 
(dashed line) in comparison to the acceleration rate (thick solid line) in the acceleration 
zone. Parameters are: L bo i,OB = 10 5 L©, er = 10 eV, Mob = lO" 6 M yr _1 , M WR = 10M O b, 
vob = 3768km/s, D = 10 14 cm, B s = 100G, x b ~ 0.2AD, B w 0.5G, n phiT w 10 n cm- 3 , 
N H ^5- 10 7 cm- 3 , n a = 1.6 • 10 20 cm 2 s- 1 , V ss 1884 km/s, T w 1127 s, r w 8.5 • 10 n cm. 
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Fig. 3. — Steady-state electron spectrum in the acceleration region for D — 5 • 10 13 cm 
(solid line), 10 14 cm (dotted line), 2 • 10 14 cm (dashed line), 5 • 10 14 cm (dashed-dotted 
line) and 10 15 cm (dashed-tripple dotted line). Parameters are: Q = lcm~ 3 s _1 , 
B w 1.3,0.5,0.30.1 and 0.05 G, n phiT w (44, 11, 2.8, 0.4 and 0.1) x 10 10 cm- 3 , N H w 
(21,5.2, 1.3,0.2 and 0.05) x 10 7 cm~ 3 , w OB = 3537, 3768, 3884, 3954 and 3977km/s, V w 
1768, 1884, 1942, 1977 and 1988 km/s, r « (9,8.5,8.2,8.1 and 8 x 10 n cm. All other pa- 
rameters are the same as used in Fig. 2. 
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Fig. 4. — Steady-state electron spectrum in the acceleration region assuming n a = = 
4 ■ 10 19 cm 2 s _1 and for D = 5 • 10 13 cm (solid line), 10 14 cm (dotted line), 2 • 10 14 cm (dashed 
line), 5 • 10 14 cm (dashed-dotted line) and 10 15 cm (dashed-tripple dotted line). All other 
parameters are the same as used in Fig. 3. 
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Fig. 5. — Energy loss rates in the convection zone due to inverse Compton scattering in the 
Thomson regime (TL, thin solid lines) and the Klein-Nishina regime (KN). The parameters 
are the same as in Fig. 2. 
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Fig. 6. — Steady-state electron spectra in the convection region at distance Ar = 10 n cm 
(dotted line), 10 12 cm (dashed line), 10 13 cm (dashed-dotted line) and 5 ■ 10 13 cm (dashed- 
triple dotted line) from the transition point r as compared to the electron spectrum in 
the acceleration region (solid line) and for a binary separation of D = 10 14 cm. All other 
parameters are the same as in Fig. 2. 
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Fig. 7. — Same as in Fig. 6 but for a binary separation of D — 10 15 cm and Ar = 10 n cm 
(dotted line), 10 12 cm (short dashed line), 10 13 cm (dashed-dotted line), 10 14 cm (dashed-triple 
dotted line) and 10 15 cm (long dashed line). 
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Fig. 8. — Steady-state proton spectra in the acceleration region for binary separations D = 
5 • 10 13 cm (solid line), 10 14 cm (dotted line), 2 • 10 14 cm (dashed line), 5 • 10 14 cm (dashed- 
dotted line) and 10 15 cm (dashed-tripple dotted line). Qq is set to unity, T e = 10 8 K. All other 
parameters are the same as in Fig. 3. 
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Fig. 9. — Steady-state proton spectra in the convection region at distance Ar = 10 11 cm 
(dotted line), 10 12 cm (dashed line) and 10 13 cm (dashed-dotted line) from the transition 
point r as compared to the proton spectrum in the acceleration region (solid line) and for 
a binary separation of D = 10 14 cm. All other parameters are the same as in Fig. 8. 
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Fig. 10. — IC spectra from the acceleration region for D = 10 14 cm and 0^ = 0° (lower solid 
line), 30° (dashed line), 60° (dotted line), 90° (dashed-dotted line), 120° (dashed-triple dotted 
line) 150° (long dashed line) and 180° (upper solid line). The normalization corresponds to 
a injected power of particles with energy E = 1 MeV prior to acceleration of ~ 0.031% of 
the total OB-wind luminosity, the emitting volume is ~ 5 • 10 37 cm 3 . For L = 30° the total 
volume-integrated (i.e. acceleration plus convection zone) IC spectrum is also shown. All 
parameters except for D are the same as used in Fig. 3. The thickness d of the emission 
region is estimated by d ~ x = x OB (Eichler & Usov 1993). 
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Fig. 11. — IC spectra from the acceleration region for D = 10 14 cm, inclination angle % = 45° 
and $ B = 0° (solid line), 90° (dotted line), and 180° (dashed line). For $ B = 180° the total 
volume-integrated (i.e. acceleration plus convection zone) IC spectrum is also shown. We 
define $ B = for the WR-star being in front of the OB-star. All other parameters are the 
same as used in Fig. 10. 
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Fig. 12. — Relativistic bremsstrahlung spectra from underlying electron distributions as 
presented in Fig. 3 for various binary separations D — 5 • 10 13 cm (dashed-triple-dotted line), 
10 14 cm (solid line), 2 • 10 14 cm (dashed-dotted line), 5 ■ 10 14 cm (dashed line) and 10 15 cm 
(dotted line), emitted from the acceleration region. The normalization corresponds to a 
injected power of particles with energy E = 1 MeV of ~0.14%, 0.031%, 0.007%, 0.001%, 
and 0.0003%, respectively, of the total OB-wind luminosity. The emitting volume varies 
between 3 • 10 37 -5 • 10 38 cm . The upper solid line corresponds to the sum of bremsstrahlung 
emission from acceleration and convection region (up to Ar = 10 13 cm) for a binary separation 
of D = 10 14 cm. 
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Fig. 13. — 7r°-decay 7-ray spectra from underlying proton distributions as presented in 
Fig. 8 for various binary separations D — 5 • 10 13 cm (dashed-triple-dotted line), 10 14 cm 
(solid line), 2 • 10 14 cm (dashed-dotted line), 5 • 10 14 cm (dashed line) and 10 15 cm (dotted 
line), emitted from the acceleration region. The normalization corresponds to a injection 
power of ~0.14%, 0.031%, 0.007%, 0.001%, and 0.0003%, respectively, of the total OB-wind 
luminosity in form of thermal particles of energy -Eo,kin — 1 MeV. The emitting volume 
varies between 3 • 10 37 — 5 • 10 38 cm 3 . The upper solid line corresponds to the sum of 7r°-decay 
photon emission from acceleration and convection region (up to Ar = 10 13 cm) for a binary 
separation of D — 10 14 cm. 
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Fig. 14. — 77 pair production opacity r 77 for a binary separation of D = 10 cm and for 
angles 9^ = 0° (lower solid line), 30° (long-dashed line), 60° (dashed-triple-dotted line), 90° 
(dashed-dotted line), 120° (dashed line), 150° (dotted line), 180° (upper solid line), and 
fixing r at 10 12 cm. 
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15. — 77 pair production opacity r 77 for a binary separation of D = 10 14 cm and varying 
10 n cm (dotted lines), r = 10 13 cm (dashed lines), and for angles 6 L = 0°, 45°, 90°, 135°, 
(from lower to upper curves). 
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Fig. 16. — 77 pair production opacity r 77 for an angle of 6^ = 45° and for a separation 
between the collision region and the OB-star surface of tob = 2.4 • 10 14 cm (dashed-dotted 
line), 1.2 • 10 14 cm (dashed line), 4.8 • 10 13 cm (dotted line), 2.4 • 10 13 cm (solid line) and 
1.2 ■ 10 13 cm (dashed-triple-dotted line) and fixing r at 10 12 cm. 
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Fig. 17. — Steady-state electron spectra for WR 140 at orbital phases 0.955, 0.2, 0.671 and 
0.8, Q = 1. See text for parameters and discussion. 
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Fig. 18. — IC spectra for WR 140 at phases 0.955, 0.2, 0.671 and 0.8 from electron spectra 
as shown in Fig. 17. The spectral changes from 7-ray absorption are not shown here. The 
EGRET 2a upper limit (Miicke & Pohl 2002) is based on observations that correspond to 
a superposition of orbital states rather determined by periastron phase, and is therefore 
compared to phase 0.955 here. 
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Fig. 19. — Broad band SED of WR 140 at orbital phase 0.955 (close to periastron). The 
injected energy in electrons and protons is equal and ~ 10~ 5 times the kinetic OB-wind 
energy. See text for further parameters and discussion. 7-ray absorption alters the spec- 
trum significantly only above ~ 60 GeV. The escaping photon spectrum remains therefore 
unchanged at this orbital phase. 
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Fig. 20. — Broad band SED of WR 140 at orbital phase 0.671. An injected electron energy 
of 3 • 10 33 erg/s in total are required to account for the synchrotron flux at this phase. For the 
protons we used QoE o kin ~ 8 ■ 10~ 7 L W (-Eo,kin = 1 MeV). See text for further parameters and 
discussion. 7-ray absorption due to photon collisions in the radiation field of the OB-star 
has been taken into account for all radiation processes: IC (solid line) and bremsstrahlung 
(dotted lines) and 7r°-decay 7-rays (dashed line). The corresponding upper and lower curves 
belong to the unabsorbed and absorbed fluxes, respectively. The cutoff in the IC spectrum 
has been extended (dotted line) using an exponential shape to guide the reader's eye. See 
text for details. 
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Fig. 21. — Steady-state electron spectrum for WR 147, Qo — 1. See text for parameters and 
discussion. 




Fig. 22. — IC spectra of WR 147 for orbital phases 0, 0.25, 0.5 and 0.75 for an underlying 
electron spectrum as shown in Fig. 21, neglecting any eccentricity of the system and assuming 
i = 90°. See text for further parameters and discussion. 77 pair production absorbes not 
more than <0.3% (> 50 GeV) and <18% (> 100 GeV) of the produced flux at orbital phases 
0.25 and 0.5, respectively (not shown in figure). No absorption takes place at phase 0. 
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Fig. 23. — Broad band SED of WR 147 at orbital phase 0. See text for further parameters 
and discussion. 




Fig. 24. — Scattering rate for e = 10 7 , 7 = 10 5 , and varying ji, p h- For curve i the scattering 
angle is n — 0.2i. 



